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Abstract 

We extend our earlier shear-transformation-zone (STZ) theory of amorphous plasticity to in- 
clude the effects of thermally assisted molecular rearrangements. This version of our theory is a 
substantial revision and generalization of conventional theories of flow in noncrystalline solids. As 
in our earlier work, it predicts a dynamic transition between jammed and flowing states at a yield 
stress. Below that yield stress, it now describes thermally assisted creep. We show that this theory 
accounts for the experimentally observed strain-rate dependence of the viscosity of metallic glasses, 
and that it also captures many of the details of the transient stress-strain behavior of those mate- 
rials during loading. In particular, it explains the apparent onset of superplasticity at sufficiently 
high stress as a transition between creep at low stresses and plastic flow near the yield stress. We 
also argue that there are internal inconsistencies in the conventional theories of these deformation 
processes, and suggest ways in which further experimentation as well as theoretical analysis may 
lead to better understanding of a broad range of nonequilibrium phenomena. 
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I. INTRODUCTION 

In the first paper of this series Q], we showed that energetic constraints determine the 
principal ingredients of a shear-transformation-zone (STZ) theory of amorphous plasticity. 
That analysis pertained strictly to the behavior of noncrystalline solids well below their 
glass temperatures. We turn our attention here to the roles played by thermal fluctuations, 
specifically, to the ways in which glassy materials make transitions from thermally activated 
creep to viscoplastic flow near yield stresses. 

The atomic mechanisms of plastic deformation are most often described as arising from 
dislocation motion. This picture breaks down in amorphous solids in which the dislocation, 
being a lattice defect, ceases to provide a useful description of the microstructural dynam- 
ics. In this paper, we describe further progress in the shear-transformation-zone theory of 
amorphous plasticity that we originally constructed to relate plastic deformation to specific 
microstructural degrees of freedom in non-crystalline solids. From its inception, our STZpic- 
ture has been based on the ideas of Cohen, Turnbull, Spaepen, Argon and others Q, El 0,0], 
who postulated that plastic deformation in amorphous materials occurs at localized sites 
usually called flow defects. A number of computational studies (e.g. jf], [?!) have provided 
support for the idea that a model based on localized defects can capture the dynamics of 
deformation in such systems. The basic premise of our version of the STZ theory is that 
these defects must be dynamic entities that carry orientational information; and our most 
striking conclusion is that, once these orientational degrees of freedom are taken into ac- 
count, the system so described exhibits an exchange of dynamic stability between jammed 
and flowing states at a stress that we identify as a yield stress. 

Earlier defect theories of deformation in glassy materials appear to us to be incomplete in 
important respects and, in some cases, to contain physically unrealistic assumptions. These 
theories generally start by assuming that the plastic strain rate e pl is the product of the 
density of defects n times an Eyring rate factor [8]: 

^ = 2n , exp (_0) sinh (_^_), (L1) 

where v is a molecular vibration frequency, AG is an activation barrier, is Boltzmann's 
constant, T is the temperature, Q is an atomic volume, and s is the deviatoric stress (i.e. 
the shear stress). The authors of these theories then attempt to describe the deformation 
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dynamics by postulating equations of motion for n. The most common choice for such an 
equation of motion has the form 



Here, k r is a thermally activated rate factor, n eq is the thermal equilibrium density of flow 
defects in the absence of external driving forces, and P is a production rate that vanishes 
when the strain rate is zero. An important example of the use of this equation is the paper 
by De Hey et al. [^]. As we shall see, our STZ theory has many features in common with 
Eq. ()1.2|) ; indeed, some form of each of the terms in this equation will appear here. 

There are important differences, however. Theories based on Eqs. and (jl.2j) make no 
attempt to describe what actually happens when a zone undergoes a shear transformation. 
Despite the fact that the Eyring formula in Eq. (jl.lj) describes the balance between forward 
and backward transitions in some kind of two-state system, the dynamics of atomic scale 
structural rearrangements described by Eq. ljl.lj) is decoupled from the population dynamics 
described by Eq. (jl.2|) . The latter equations, therefore, implicitly assume that there exists 
some fast relaxation mechanism - faster than any other rate introduced explicitly in the 
theory - which causes zones instantaneously to lose their memory of prior transformations. 

Another major difference is in the choice of the production function P. In some theories, 
P is chosen to be linearly proportional to e pl , which is impossible because P must be a 
non-negative scalar while e pl is a tensor that can change sign. Such theories generally are 
used only in cases where e pl is positive, which may mean that the authors intend to use 
the magnitude of e pl in more general situations; but the latter convention also would be 
unsatisfactory because \e pl \ is a nonanalytic function that is not likely to arise from any 
first-principles analysis of molecular mechanisms. In our own earlier paper Q], we tried 
using the rate of plastic work, se pl , in our analog of the production term. That function is 
a scalar with a satisfactory physical interpretation but, as we remarked there, it also suffers 
from a sign problem because it can be negative during unloading. We believe that we have 
solved this problem in 1], and shall make extensive use of the technique described there in 
what follows. 

Spaepen DO has introduced an important modification of the above ideas by postulating 
that the defect density n is directly determined by the excess free volume in the system, v /, 
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via a relation of the form 

( v *\ 

n oc exp , (1.3) 

where V* is a molecular volume. Spaepen's proposal is that the production term P in 
Eq. (jl.2j) should be proportional to the growth rate of Vf. More recently, Johnson et al.jlj| 
have proposed a dynamic free volume model that includes a phenomenological parameter 
very roughly analogous to the yield stress that emerges from our STZ theory. Unfortunately, 

n 

Spaepen, Johnson and others (e.g. [9|) postulate an equation of the form Vf oc e pl , which 
again violates symmetry or analyticity requirements. In short, we believe that these theories 
require a critical reformulation in order to enable a meaningful atomic-scale analysis of 
amorphous plasticity. 

In addition to reformulating earlier theories of amorphous plasticity, one of our principal 
goals in this paper is to gain as simple as possible an understanding of recent experimental 
results on plastic flow in metallic glasses. In particular, we refer to work by Kato et al.fl^ 
on amorphous Pd 40 Ni 10 Cu 30 P2tb an d the results of Lu et al.[13], who measured properties 
of bulk amorphous Zr 4 i. 2 Tii 3 8 CU12.5 Nii Be22.5 over a remarkably wide range of strain rates 
and temperatures. 

For the sake of simplicity, we restrict ourselves throughout this paper to the quasilinear 
theory discussed in [l|. We also make other assumptions which, as we shall indicate at 
several places in our presentation, raise fundamental issues regarding nonequilibrium states 
in deforming solids. 

Our order of presentation is as follows. In Section II, we review the elements of the 
low-temperature, quasilinear theory [l] in a way that prepares concepts and notation for 
the discussion of thermal relaxation in Section III. Section IV is devoted to analysis of 
these theoretical results and comparison with experiments. Finally, in Section V, we return 
to various fundamental issues that are raised earlier in the paper or are relevant to its 
conclusions. In particular, we discuss how our original, fully nonlinear STZ theory may 
need to be invoked in order to make better contact with atomic-scale mechanisms; we briefly 
address the important issue of shear banding; and we discuss the question of departures from 
thermodynamic equilibria in driven systems - a question that we shall not be able to avoid 
in the present analysis. 
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II. ELEMENTS OF THE LOW-TEMPERATURE STZ THEORY 



In order that this account be reasonably self-contained, and to provide some new per- 
spectives about the discussion that follows, we start by reviewing the basic elements of the 
low-temperature STZ theory. 



Throughout the following analysis, we take our original STZ picture 10| more literally 
than perhaps is necessary. We assume that, instead of being structureless objects as in 
the earlier theories described in the Introduction, the STZ's are two-state systems. In 
the presence of a shear stress, they can deform by a finite amount in one direction before 
becoming jammed and, when jammed in one direction, they can transform in the opposite 
direction in response to a reversed stress. Importantly, our STZ's are ephemeral; they are 
created and annihilated during irreversible deformations of the material. 

The literal interpretation to be used here requires that all STZ's have approximately 
the same size and dynamic properties. To visualize such an STZ, think of a void in an 
elastic material, and place a small group of molecules inside it in such a way that their 
average free volume is somewhat larger than that for most other molecules in the system. 
The void has some degree of structural stability; it can deform elastically but, because 
of the configuration of molecules on its surface, it resists collapse. Rearrangements of the 
molecules that are caged within the void couple to its shape and, therefore, to the stress field 
in the elastic medium in which the void is embedded. This picture suggests that the system 
undergoes two distinct kinds of irreversible events: volume conserving shear deformations, 
i.e. the STZ transformations, and dilations or contractions in which the STZ's are created or 
annihilated. (Nothing in this picture rules out the possibility that transient dilations occur 
during intermediate stages of the shear transformations.) 

For simplicity, as in ll. Ilpj. we consider a two-dimensional model, and we subject this 
system only to pure shear deformations. As we shall show at the beginning of Section 
IV, the properties of this model are easily reinterpreted in terms appropriate to simple 
uniaxial stress experiments in three dimensions, so long as we are willing to assume that 
the system remains spatially uniform. We further restrict ourselves to molecular materials 
in contact with thermal reservoirs (as opposed to granular materials or foams), so that 
we may assume that an ambient temperature determines an underlying fluctuation rate. 
At low temperatures, thermal fluctuations provide an attempt frequency for stress-induced 
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molecular rearrangements, but are too small to activate transitions over energy barriers in 
the absence of external driving forces. This attempt frequency will never appear explicitly 
in our analysis but, rather, is embedded in the rate factors to be introduced below. 

For present purposes, we need to consider only situations in which the orientation of the 
principal axes of the stress and strain tensors remains fixed. That is, we do not consider 
situations in which a fully tensorial version of the STZ theory will be necessary, as in the 
necking calculations reported in Q]. Therefore, it is sufficient to assume that the population 
of STZ's consists simply of zones oriented along the two relevant principal axes of the stress 
tensor. Exactly the same equations as the ones we shall use here can be derived starting 
from the assumption that the a priori distribution of orientations of the zones is continuous 
and Symmetri e.R0 Without .oss of geaeraUt, therefore, we let the deviate stress be 
diagonal along the x, y axes; specifically, let s xx = —s yy = s and s xy = 0. Then choose the 
"+" zones to be oriented (elongated) along the x axis, and the "— " zones along the y axis; 
and denote the population density of zones oriented in the "+/— " directions by the symbol 
n±. 

With these conventions, the plastic strain rate is: 



Here A is a material-specific parameter with the dimensions of volume (or area in strictly 
two-dimensional models), which must be roughly the same order of magnitude as the volume 
of an STZ, that is, a few cubic or square atomic spacings. The remaining factor on the right- 
hand side of Eq. (|2.1|) is the net rate per unit volume at which STZ's transform from "— " 
to "+" orientations. R{s)/r and R(-s)/t are the rates for "+" to "-" and "-" to "+" 
transitions respectively, where r is the time scale that characterizes the low-temperature 
plastic response. For simplicity, we write these rates as explicit functions of only the devia- 
toric stress s, although they depend implicitly on the temperature and pressure and perhaps 
other quantities. 

The equation of motion for the populations n± generally must be a master equation of 
the form: 



The first two terms on the right-hand side are the stress-driven transition rates introduced 
in the preceding paragraph. There is no analog of these terms in Eq. ([1.2)1 . They de- 




(2.1) 




(2.2) 
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scribe volume-conserving, pure-shear deformations which preserve the total population of 
the STZ's. The last two terms in parentheses, proportional to T, describe creation and 
annihilation of STZ's. In the low-temperature theory, T is nonzero only when the plastic 
strain rate is nonzero; the molecular rearrangements required for creating or annihilating 
STZ's cannot occur spontaneously, that is, in the absence of external driving forces 

The assumption in Eq. ()2.2|) that the annihilation and creation rates are both proportional 
to the same function T has serious implications in this theory. Among those implications is 
the requirement that strain-rate independent constant. Note that is the total 

density of zones generated in a system that is undergoing steady plastic deformation. It is 
not the same as the quantity n eq introduced in Eq. (jl.2|) . which is the equilibrium density at 
nonzero temperature and zero strain rate, and ordinarily is said to go rapidly to zero as the 
temperature decreases below the glass transition. On the other hand, property of 

low-temperature materials at non-zero strain rates. 

The form in which we have cast Eq. (j2.2|) is consistent with a fundamental assumption 
that we are making about the nature of our low-temperature theory. Specifically, we are 
assuming that the only relevant time scales at low temperatures are To and the inverse of the 
strain rate. This means that, under steady-state conditions at strain rates less than some 
value of order t _1 , the number of events in which the molecules rearrange themselves is not 
proportional to the time but to the strain. That picture seems intuitively reasonable. If 
the system requires a certain number of STZ-like rearrangements in order to achieve some 
deformation, then it should not matter (within limits) how fast that deformation takes place. 
The picture breaks down, of course, when there are competing rearrangement mechanisms. 
For example, we shall see that the density of STZ's becomes strain-rate dependent when 
we introduce thermal fluctuations, because such fluctuations will induce rearrangements at 
a rate that is independent of the strain rate. We also expect that the picture may fail in 
polymeric glasses or polycrystalline solids, where more complex components may introduce 
extra length and time scales. 

This simple dimensional argument, leading to a nonzero, rate-independent value of n^, 
already hints at a fundamental difficulty in theories of the kind summarized by Eqs. (jl.l|) 
and ()1.2|) . These theories have no sensible low-temperature limit because both the strain 
rate in Eq. (jl.l|) and the rate factor k r in Eq. (jl.2|) vanish rapidly as T — > 0. Yet, even at 
temperatures so low that thermal fluctuations cannot cause molecular rearrangements, such 
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systems must deform irreversibly when sheared. 

We shall use the energetic arguments introduced in to determine the factor T in 
Eq. (|2.2|) . but first we must discuss the state variables and specific forms for the transition 
rates. We define dimensionless state variables by writing 

A= n+ + n -, A= n +~ n -. (2.3) 



In a more general treatment jl^ . Ilf! Il7 |. A remains a scalar density, but A becomes a 
traceless symmetric tensor with the same transformation properties as the deviatoric stress. 
We also define: 

S = ~(R(-s)-R(+s)), C = ±(r(-s)+R(+s)), T=~. (2.4) 
Then the STZ equations of motion become: 

r e pl = e C(s)(AT(s)-A); (2.5) 



r A = 2 C(s) (A T(s) - A) - r(s, A, A) A; (2.6) 

and 

r A = r(a,A,A)(l-A). (2.7) 

Here, eo = Xn^ is roughly the fraction of the total volume of the low-temperature system 
in steady-state flow that is covered by the STZ's. This is a material-specific quantity. If eo 
is small, then the disorder induced in the system by deformation is small. Conversely, if eo 
is large, then the STZ-like defects cover the system and the material in some sense "melts" 
under persistent straining. 

Throughout this paper, we shall use only what we call the "quasilinear" version of these 
equations. 3] That is, we note that T(s) and C(s) are, respectively, antisymmetric and 
symmetric dimensionless functions of s, and write: 

T( 8 ) = - = s; C(s) = 1, (2.8) 

Sy 



where s y will turn out to be the yield stress. The choice C, „ , 



s) = 1 is, in effect, our definition 
of the time constant To. As pointed out in earlier papers llSf, this quasilinear approxima- 
tion has important shortcomings. Neglecting the stress dependence of C(s) means that we 
overestimate the amount of plastic deformation that occurs at small stresses and therefore 
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also overestimate the rate at which orientational memory disappears in unloaded systems. 
Moreover, the quasilinear approximation is too simplistic to be related directly to atomic 
mechanisms, a point that we shall comment upon further in Section V. Nevertheless, the 
quasilinear theory has the great advantage that it is mathematically tractable and easy to 
interpret. It will serve to illustrate the main points that we wish to make in this paper, but 
aspects of the nonlinearities associated with C and T will need to be reintroduced before we 
shall be able to understand fully the nonequilibrium behavior of amorphous solids. 
Eqs. (|2.5N2.7jl now become: 

T e pi = e (As~-A); (2.9) 
t A = 2 (A 5- A) -T(s,A,A) A; (2.10) 

and 

r A = r(s,A,A)(l-A). (2.11) 

The quantity T/t is the STZ creation rate and therefore plays exactly the same role 
as the phenomenological defect production rates that we discussed - and complained about 
- in the Introduction. We can derive an expression for that rate by using the energy-balance 
argument introduced in [lj. As before, we start by writing the first law of thermodynamics 
in the form: 

22*8= ( A S - A) S = e s y ^ ip(A, A) + Q{§, A, A). (2.12) 

The left-hand side of Eq. (|2.12|) is the rate at which plastic work is done by the applied stress 
s = s y s. On the right-hand side, eo s y ip is the state-dependent recoverable internal energy, 
and Q is the dissipation rate. So long as the STZ's remain uncoupled from the heat bath, Q 
must be positive in order for the system to satisfy the second law of thermodynamics, that 
is, for the work done in going around a closed cycle in the space of variables s, A, and A to 
be non-negative. 

As argued in [l|, the simplest and most natural choice for T - and, so far as we can tell, 
the only one that produces a sensible theory - is that it be the energy dissipation rate per 
STZ. That is, 

Q(S,A,A) = ^Ar(S,A,A), (2.13) 

TO 

With this hypothesis, we can use Eqs. ()2.10|) and (J2.11|) to write Eq. (j2.12|) in the form 

2(AS-A)S=^r(l-A) + ^(2(Aa-A)-rA)+Ar. (2.14) 
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Then, solving for T, we find: 

r = 



2 (As — A) (S-dip/dA) 



A + (1 - A)(&ip/dA) - A (d^ /OA) ' 
To assure that T remains non-negative for all s, we must let 



(2.15) 



di) _ A 
dA ~ A ! 



(2.16) 



so that the numerator becomes 2 A (3 — A/A) 2 . Then (see |l[), we choose 



0: 



so that 



^(A,A) 



r(S,A,A) 



A 



1 + 



A 5 



a 2 ; ' 

4 A (AS - A) 2 



(2.17) 



(2.18) 



1 + A) (A 2 - A 2 )' 

This result has the physically appealing feature that it diverges when A 2 approaches its 
upper limit A 2 , thus enforcing a natural boundary for dynamical trajectories in the space of 
the state variables A and A. 

It is convenient at this point to replace the variable A by m = A/A, so that the equations 
of motion become: 



T e pl = e A(s-m); 



Tq fa = 2 (s — m) 



1 - 



2 m (s — m) 
(l+A)(l-m 2 ) 



and 



A 

To A 



4(s — m) 
1 — m 2 



1- A 
1 + A 



At the stable fixed point of Eq.flEZU, A = 1, Eq.flODJ) becomes 

2 (s — m) (1 — sm) 



r m 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



(1-m 2 ) 

which exhibits explicitly the exchange of stability at s — 1 between jammed states with 
m = s and flowing states with m = 1/5. 



III. THERMAL EFFECTS 

We return now to Eq. (J2.2j) . the low-temperature master equation for the STZ population 
densities n±, and ask what changes need to be made in order to incorporate thermal effects 
at temperatures near the glass transition. One obvious possibility is to modify the rate 
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factors R(±s) to include thermal activation across energy barriers; indeed, we eventually 
shall have to do that. (See Section V.) However, our quasilinear approximation makes it 
difficult to do this systematically. 

The more important thermal effects are those that are completely missing in Eq. (|2.2|) . 
specifically, the thermally assisted relaxation - i.e. aging - of the STZ variables that can 
occur spontaneously in the absence of external driving or plastic strain rate. There are two 
ways in which relaxation must occur. First, thermal fluctuations ought to act much like 
deformation-induced disorder in causing the n± to relax toward their steady-state values 
rioo/2. Second, there should be some annealing mechanism that causes the total STZ popu- 
lation to decrease. Both of these mechanisms involve dilations and contractions of the kind 
associated with creation and annihilation of STZ's; thus, again for simplicity, we assume 
that there is just a single relaxation rate, denoted p(T)/tq, that characterizes them. As we 
shall see, that rate will have the Vogel-Fulcher or Cohen-GrestQ form, rapidly becoming 
extremely small as the temperature T falls below the glass temperature. Specifically we 
expect that p{T) has the form 

/ Ay dil \ 

p{T) = p exp l-—^ 1 , (3.1) 

where po is a dimensionless prefactor, AV dd is the activation volume required to nucleate 
a dilational rearrangement, and v/(T) is usually identified as the free volume. The Cohen- 
Grest approximation for Vf(T) has the form 

where v , T and Ti are fitting parameters. This expression was found by Masuhr et al. 
(^(J and Lu et al. [kJ to provide a fairly accurate fit to their data. (We shall not use this 
formula explicitly in what follows.) 

In accord with the preceding remarks, our proposed form for the modified master equation 

is: 

r h± = R(^s)n T — R(±s)n± 

+ (r( S , A, A) + p(T)) {^f - n ± ) - k p{T) n±. (3.3) 

The first and second appearances of p(T) on the right-hand side of Eq. (j3.3J) correspond, 
respectively, to its two roles described above. The second of these terms, the quadratic form 
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r-T + V(T-T ) 2 + Tir, (3.2) 



with a dimensionless multiplicative constant k, is similar to the n 2 term on the right-hand 
side of Ea. (jl.2j) . This bimolecular mechanism has been discussed extensively in references 
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231. 



Equations (j2.9H2.llf) now become: 



r„ 6* = e„ (As - A) 



(3.4) 



(unchanged from before); 



t A = 2 (As- A) - (r(s, A, A) + p(T) ) A - K p(T)AA; 



(3.5) 



and 

r A = (T(s, A, A) + p(T)) (1 - A) - k p(T) A 2 . (3.6) 

The next step is to repeat the energy-balance analysis of Eqs. (12.121 - 12. 18|) to recompute 
the function r(s, A, A). We assert that Eq. ()2.13|) relating T and the dissipation rate Q must 
remain unchanged by the addition of the thermal relaxation terms in Eq. (l3.3|) . That is, V 
in Eq. (j3.3|) must be the energy dissipated per STZ when plastic work is done on the system. 
The expression for the internal energy ip(A,A) must remain as given by Eq. (J2.17|) because 
Eq. (|2.16|) is still required by the non-negativity condition. The result, after inserting the 
terms proportional to p(T) into Eqs. ()2.14|) and ()2.15|) . and transforming to m = A/A, is: 

"4 (s - m) 2 + 2 p(T) + k p{T) A (1 + m 2 ) 



r(S,A,m) + p(T) = A 



;i + A)(l 



Af(5,A,m,T). 



(3.7) 



Note that V is non-negative, as is necessary because it is a prefactor for the annihilation and 
creation rates. The non-negativity condition no longer strictly applies to V itself because 
the system is now coupled to the heat bath. 
The new equations of motion for m and A are: 



and 



T e pl = e A (s - m) 



r rh = 2 (s — m) — m T(s, A, m, T), 



r j = f(S, A, m, T) (1 - A) - « p(T) A. 



(3.8) 
(3.9) 

(3.10) 
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As in previous presentations, we recover Bingham-like plasticity for small p and for 
stresses above s y , i.e. s > 1. In that case, the equations of motion revert to Eqs. (|2.19l 
I2.21|) so that, in steady state, A — > 1 and m — > 1/5. Thus, 



■pi . e o /- 1 



(?*«- 3-^ ; 5>1. (3.11) 

This theory exhibits no power-law rheology for flow above the yield stress. As we shall see, 
however, it does exhibit behavior that looks like superplasticity. 

The more interesting behavior is thermally assisted creep below the yield stress, s < 1, 
and the transition between creep and flow near s = 1. For present purposes, we only need to 
consider cases in which p(T) <C 1. In the steady-state creep region where s<Cl and m = s, 
Eqs. fl3.7|) and ()3.9|) tell us that s — m is small of order mp. Then the quantity (3 — m) 2 in 
the numerator of T is negligible compared to the other terms. The steady-state version of 
fl3.1()|) becomes simply 

1 - A n (k) « kA 2 n {k); A n {k) = — (v / TT4^ - l) , (3.12) 

where the subscript N denotes the low-stress Newtonian limit. 

In the complete absence of a driving stress 3, the exact steady-state solutions of Eqs. (J3.9j) 
and ()3.1()|) are m = and A = Ajy(/c). Thus, the analog of n eq in Eq. (|1.2jl is 11^ An(k), which 
is a temperature-independent quantity instead of being, as is usually assumed, a rapidly 
varying function of the form exp(— AV dd /vf(T)). Indeed, the temperature dependence of 
is often invoked in the context of equations such as Eq. (jl.2j) to interpret calorimetric 



n. 



<-■<-} 



data 



24j Because of its lack of a temperature-dependent n^, our present theory cannot 
predict the specific heat peak near the glass transition that is seen by the latter authors. We 
could fix this problem in an ad hoc manner by assigning some temperature dependence to 
and/or k. However, such a procedure would simply gloss over the fundamental difficulty that 
we are facing here - that the limiting steady-state value of A must depend upon the order 
in which we take the limits T — ► and e pl —>■ 0, but that no such behavior appears in our 
equations. This is the same situation that we discussed in the paragraphs following our first 
introduction of in Eq. (j2.2|) : we shall return to it in Section V. For the present, we leave 
the situation as is, with a temperature-independent Ajv, and with the understanding that 
we cannot yet use these equations to describe behavior much above the glass temperature. 
To compute the viscosity, it is easiest first to set m = A = in Eqs. (j3.9|) and (|3.10j) and 
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FIG. 1: Theoretical curves of tensile stress versus strain for the bulk metallic glass 
Zr4i.2 Tii3.8 Cui2.5 Niio Be22.5 at several different temperatures as shown. The strain rate is 
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For clarity, the curves have been displaced by constant increments along 



the strain axis. 



eliminate T to find 



(s-m) (1 - A) 



K P (T) 



m A, 



(3.13) 



which is an exact relationship between the steady-state values of s, A and m at any strain 
rate. Combining these results in the limit s ~ m — ■> 0, we compute the Newtonian viscosity 

VN = lim -i- = (3.14) 
which confirms our expectation that p(T) is the rate function that governs viscous relaxation. 



IV. ANALYSIS AND COMPARISON WITH EXPERIMENTS 

Before making comparisons with experiments, we must return to the question of how 
to generalize our two-dimensional theory into one that can be applied to three-dimensional 
experiments. The basic structure of our equations of motion must be preserved, with due 
respect for the relevant symmetries, in any generalization of this theory to higher dimensions. 
That is, our variables s, e pl and m must become traceless symmetric tensors, and A must 
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remain a scalar, so that Eqs. ()3.8j) and ()3.9|) become: 

r eg = e A (% - my), (4.1) 

and 

r my = 2 (sy - m i3 -) - f(s, A, m, T) my. (4.2) 

The energy-balance analysis yields 

rfz a TO T ^ - 2 ^ ~ ~ + 2 ^( T ) + ^( T ) A ( x + 

(1 + A) (1 — m z ) 

where we are using the summation convention, and m 2 = (1/2) my my. 

The experiments in which we are interested involve only uniaxial stresses, say, in the x 
direction. If there are no tractions in the y or z directions, and if we can ignore spatial 
nonuniformities, then the total stress cry is diagonal with a xx = a, a yy = o zz = 0. Similarly, 
for the deviatoric stress, s yy = s zz = —s xx /2; and m yy = m zz = —m xx /2. Thus, m 2 = 
3/4 m xx . We can now make the xx components of Eqs. (|4.1|) - (|4.3J1 look exactly like Eqs. 
(I3~£l) . (jHU), and flHHI) by defining m 2 = m 2 , m = ^3/4m xx , s = yfe/Zs^, e pl = e p J x , 



and by replacing eo by e' = y4/3eo- Note that the dynamical exchange of stability, i.e. 
plastic yielding, still occurs at s = 1; thus our scaling of the stress by s y remains correct. 
When comparing to the data, which is presented in terms of the uniaxial stress a = a xx , we 
write a = (3/2)s xx = \f3s y s or, equivalently, a = a y s, where a y is the tensile yield stress. 
Equation ()3.14|) . however, remains unchanged with e rather than e' , because viscosity now 
is expressed as rj = s xx /2e p x l x = s/^/3e pl . 

To illustrate the principal results of our analysis, we first follow the lead of Kato et al 
jl^ and Lu et al by looking for scaling in the steady-state behavior of our system. 
To be specific about what we mean here, we show in Figs. 1 and 2 sets of stress-strain 
curves for various temperatures and fixed strain rates. As we shall explain shortly, these 
figures are to be compared with Figs. 1 and 2 in (See also Fig. 1 in and Fig. 1 



m ll2j.) A general feature of these curves is that, when the strain rate is held constant, 
the stress rises through a maximum, decreases as the material softens, and then reaches a 
steady-state value. We shall discuss the initial transients later in this Section, but look first 
at the late-stage, steady-state behavior. 

We compute the steady-state flow stress as a function of the strain rate by solving 

n n 

Eqs. ()3.9p and ()3.10|) with m = A = 0. Then, as in [12J and [13j], we plot s = s/s y as & 
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FIG. 2: Theoretical curves of tensile stress versus strain for the bulk metallic glass 
Zr4i.2 Tii3.g Cui2.5 Niio Be22.5 at several different strain rates as shown. The temperature is 
T = 643 K. For clarity, all but the first of these curves have been displaced by the same amount 
along the strain axis. 

function of rj^e pl for eight different values of the relaxation rate p(T) corresponding to the 
eight different temperatures for which data are reported in [13]. The results are shown in 
Fig. 3. As discovered by Kato et al [l^ . all of these curves lie on top of one another for 
stresses s < 1 but, in our case, they diverge from each other in the flowing regime, s > 1, 
where the Bingham-like behavior shown in Eq. ()3.11|) sets in. 

Figure (4a) contains the same theoretical curves as those shown in Fig. 3, but plotted 
there as tensile stress versus scaled strain rate, and compared with experimental data taken 
from Fig. 9a of |13J]. The same theoretical functions and data points are replotted in Fig. 
4b to show the normalized viscosity, r]/i]N as a function of the scaled strain rate. The latter 
figure is directly comparable to [13j, Fig. 9b. Note that the range of strain rates shown in 
Figs. 4 corresponds to the range of the experimental data and is substantially smaller than 
that shown in Fig. 3. The theoretical curve that lies above the rest at high strain rates is 
for T = 683 K, the highest of the temperatures reported in |l3j]. The data points at that 
temperature all lie at scaled strain rates that are too small to test this predicted breakdown 
of the scaling law. 

Our Fig. 5a shows individual theoretical and experimental curves of tensile stress as a 
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FIG. 3: Scaling behavior in the STZ theory: shear stress s as a function of strain rate scaled by 
rjjsf. This graph is plotted for the same set of temperatures as shown in Fig. but for a larger 
range of strain rates. 

function of (unsealed) strain rate for different temperatures. Here, the experimental data is 
from jl^, Fig. 7. These curves are replotted in Fig. 5b to show (unsealed) viscosity as a 
function of strain rate, analogous to |l3j], Fig. 8. 

In constructing these figures, we have determined our theoretical parameters as follows: 



We have used the value of the room-temperature tensile yield stress reported in 



[13], 



<t v = 1.9 GPa. Thus a = 1.9 s GPa. 



Rather than using the Cohen-Grest formula with parameters from |20l|, we have taken 

n 

the limiting (vanishing strain rate) Newtonian viscosities directly from [13J], Fig. 10, 
and have checked that these values (apart from one apparently misplaced point) are 
consistent with the data points in 13J, Fig. 9. 



Given the above constraints, we are left with only two parameters, k and €q/tq (or , 
equivalently, e' ), that can be adjusted to fit the steady-state experimental data in (l3j, 
Fig. 7. Because we know experimental values for the limiting viscosities t]n{T) at 
different temperatures, a value for the ratio ^q/tq in Eq. (j3.14j) determines the overall 
scale factor for the function p(T). Thus our fitting procedure has been to start by 
choosing values of n and e /r . We then use those values to determine p{T) and to 
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Newtonian viscosity * strain rate (Pa) 



FIG. 4: Tensile stress and viscosity as functions of scaled strain rate rmer . The data points for 

n 

the bulk metallic glass Zr4i.2 T113.8 CU12.5 Niio Be22.5 are taken from |13lj. Figs. 9a and 9b. The 
solid curves are theoretical results computed for the same set of temperatures as shown. 

compute steady-state solutions of Eqs. (j3.9j) . ()3.1U|) and (j2.19|) . From these solutions 
we compute steady-state stress versus strain-rate relations that can be compared with 
the experimental data. We then iterate this process to find best-fit values for k and 

eo/V 

Our best-fit parameter values, obtained by the procedure outlined above, are k = 120 
and eo/r = 260 s~ x . Our corresponding values of p(T), along with our estimates of the 
viscosities, are shown in Table 1. We emphasize that these values are not much better than 
order-of-magnitude estimates. 

As is obvious in Fig. 5, there is scatter in the experimental data, and there may also be 
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FIG. 5: Tensile stress (a) and viscosity (b) as functions of strain rate for different temperatures as 
shown. The data points are for the bulk metallic glass Zr4i.2 Tii3.g CU12.5 Niio Be22.5 as reported 
by Figs. 7 and 8. The solid lines are theoretical curves. 

systematic errors. For example, the two points at the highest strain rates for T = 663 K fall 
well below our theoretical curve for that temperature while the theoretical fits look good 
for the temperatures on either side, i.e. T = 643 K and 683 K. We could improve the fits to 
all three of these curves by choosing a substantially larger value of k; but we would do this 
at the expense of poorer fits at lower temperatures. However, the values of the viscosities 
that we can deduce from [3], Fig. 10 seem uncertain, possibly by factors of 2 or 3; so our 
estimates of p(T) and, therefore, all our theoretical curves - especially those at the lower 
temperatures - might be modified by more accurate viscosity data. 

Within the above uncertainties, our theoretical fits to the experimental data are relatively 
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Temperature, K 


573 


593 


603 


613 


623 


643 


663 


683 


Viscosity, PaS 


4.00 x 10 14 


4.03 x 10 13 


8.99 x 10 12 


4.03 x 10 12 


7.29 x 10 11 


4.27 x 10 10 


3.68 x 10 9 


3.99 x 10 s 


P 


1.07 x 10~ 8 


1.06 x 10- 7 


4.77 x 10~ 7 


1.06 x lO" 6 


5.88 x 10- 6 


1.01 x 10~ 4 


1.17 x 10~ 3 


1.15 x 10~ 2 



TABLE I: Experimental data for viscosity taken from 
calculations. 



131 ] . and values of p used in the present 



insensitive to our choices of the two parameters that we have allowed ourselves. On the one 
hand, this insensitivity gives us confidence in the basic structure of our theory; on the other 
hand, it means that we cannot yet test the theory in as much detail as we would like. For 
example, the yield stress s y probably ought to be a function of temperature, decreasing 
slowly (in contrast to the rapidly varying viscosity) from room temperature through the 
experimental range. Its behavior should be roughly like that of the shear modulus, which 
must soften as the system approaches the glass transition. Also, as we have discussed above, 
we expect that eo ultimately will be temperature dependent because it is proportional to 
rioo. We could improve the fits, for example, at the lower temperatures shown in Fig. 5, 
by incorporating such temperature dependences into our equations; but it seems to make 
little sense to do so without first resolving various uncertainties in both the theory and the 
experiments. 

Our main conclusion from this steady-state analysis is that we are observing a transition 
from thermally assisted creep to viscoplastic flow in the neighborhood of the dynamic yield 
stress. At low stresses and strain rates, the linear response relation contains only the factor 
t]n oc l/p(T), thus we obtain the simple scaling. Near the yield stress, however, our theoret- 
ical strain rate increases by several orders of magnitude for small increments of stress, and 
the experimental behavior tracks this trend accurately. This behavior resembles superplas- 
ticity. Interestingly, the theoretical scaling persists through the "superplastic" region and 
does not break down until true viscoplastic flow begins. 

Before returning to the transient stress-strain curves shown in Figs. 1 and 2, we show 
one other steady-state prediction of our theory for which there is no data in [13[ , but which 
seems to be potentially important. In Fig. we have used our steady-state solutions for 
m and A to plot the dimensionless internal energy ip = (A/2) (1 + m 2 ) as a function of the 
strain rate e pl , for the set of temperatures used in the preceding figures. The energy ip might 
be measured by differential scanning calorimetry, as has been done by Hasan and Boyce in 
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Strain rate (s ) 

FIG. 6: Dimensionless energy ifi as a function of strain rate for the same set of temperatures shown 
in the preceding figures. 



studies of polymeric glasses. 



25. 



De Hey et al |9| also report DSC results, which they 
interpret as measurements of free volume. We assume that Fig. 4 in |9|, apart from the 
scale on the vertical axis, is at least qualitatively the same as a graph of ip versus strain 
rate analogous to our Fig. 6. The two figures are similar to one another if we look at our 
theoretical curves only at small strain rates. 

An important feature of our Fig. 6 is that, at fixed strain rate, ip decreases as T in- 
creases. That trend is exactly what we expect for thermally assisted creep; the higher the 
temperature, the fewer STZ's are needed in order to sustain a given flow rate, and thus the 
smaller the internal energy. Note, however, that our theory predicts that all of these curves 
converge to a single, temperature independent, value i/j — *> Ajy(/c)/2 in the limit e pl — > 0. On 
the other hand, De Hey et al. plausibly assert that the equilibrium value of the STZ density, 



n. 



eq ■ 



should be an increasing function of temperature. If so, these curves must cross each 



other and the trend in temperature dependence must be reversed at small enough values of 
e pl . An extension of these measurements to smaller strain rates might therefore provide a 
test of our differing assumptions about the STZ density n eq . 

So far, we have examined only steady-state behavior. We turn next to stress-strain curves 
obtained in constant strain-rate experiments such as those shown in our Figs. 1 and 2 and 
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0. 



in Figs. 1 and 2 of .13]. To plot these curves, we solve 

¥* ^-l/f I *<•"-*">■ (4 ' 4) 

along with Eqs. (j3.9j) and (j3.10|) to com put e s as a function of the total strain. We use 
the value of Young's modulus given in E = 96 GPa, to estimate E/a y = 50. For 
time dependent calculations, we must choose a value of the time scale r , which multiplies 
the time derivatives in Eqs. (j3.9j) and (J3.10j) . To do this, we keep the ratio €q/t fixed at 
its value obtained from the steady-state calculations, and we adjust the value of eo so as 
to fit the transient stress-strain curves. Our best-fit value is eo = 0.7, which means that 
t = 2.7 x 10~ 3 s. We emphasize that this value, like our estimates for k and €q/t , remains 
highly uncertain. It is possible that we may eventually be able to sharpen our estimate of 
the time scale r by using the time- dependent stress-relaxation data shown in jl^ . In fact, 
we can reproduce those results about as well as the other results shown here, but our theory 
in this case is especially sensitive to our quasilinear approximation plus other uncertainties 
regarding the experiments. 

Our Figs. 1 and 2 are drawn so as to be directly comparable to Figs. 1 and 2 in jl^ : 
that is, we use the same strain rates and temperatures. The one other parameter that we 
must choose for solving Eqs. ()3.9j) . (j3.10j) . and ()4.4j) is the initial value of A, which we 
denote by A . Theoretically, the smallest value of A that can be achieved by annealing is 
Aat(k) = 0.087 (for k = 120); therefore we have used this value of A for these two figures. 
The initial values of m and s are always chosen to be zero. 

In all cases, the agreement between theory and experiments seems satisfactory given the 
various uncertainties. The peak heights and positions for fixed strain rate e total = 0.1 s~ l 
and varying temperatures in Fig. 1, and for fixed temperature T = 643 K and varying 
strain rates in Fig. 2, are within about ten percent of their experimental values. The 
experimental curves for low temperatures and large strain rates end where the samples 
break; the dashed lines in our figures indicate our theoretical extensions of those parts of 
the curves for which no experimental data is available. The one systematic discrepancy 
is that our initial theoretical slopes are smaller than the experimental ones. This is an 
artifact of our quasilinear approximation, which ignores the strong stress and temperature 
dependence of the factor C{s) in Eqs. (J2.5|2.6J) . (See remarks in Section V.) In the limit 
of low temperatures and very small strain rates, our theory predicts this initial slope to be 
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FIG. 7: Tensile stress as a function of strain for several different values of Ao. Curves are plotted 
for e total = 3.2 x 1(T 2 s' 1 at T = 643K. 



E/(l + E e' A(]/2a y ) instead of simply E. This does not happen in the fully nonlinear theory 

presented in [lOj. We have chosen not to include these nonlinear effects in our calculations 

here because they introduce additional undetermined parameters, and they are not necessary 

to describe the shear softening and shear thinning observed at higher stresses near the glass 

transition temperature. 

' finally, in Fig. we use the material parameters deduced above for the system studied 

in 3] to plot stress-strain curves for different Ao's, all at temperature T = and 

e total = 3.2 x 10 -2 s -1 . The different Ao's correspond to different initial states of disorder 

produced by varying the annealing times and temperatures. Presumably, annealing for 

longer times at lower temperatures produces smaller values of A ; but it seems difficult to 

make quantitative estimates of this effect. These curves may be compared qualitatively 
rn 

with those shown in |9j, Fig. 9, where larger initial densities of STZ's produce larger plastic 
responses and correspondingly smaller overshoots during the early stages of deformation. 



V. CONCLUDING REMARKS 

The comparison with experiments discussed in the preceding section leads us to believe 
that the STZ theory captures the main features of the experimental data, but that we shall 
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have to improve it in specific respects if we are to develop it into a yet more quantitative, 
predictive description of plastic deformation in amorphous solids. We conclude this paper 
by identifying three directions for the next phases of these investigations. 
Fully nonlinear, temperature dependent transition rates: When we examine the quasilinear 
STZ theory in the context of a theory that includes thermal fluctuations, we see that it is 
a special case in which the shear rearrangements are not being modeled as realistically as 
the dilations or contractions. To see this in more detail, go back to our original^], fully 
nonlinear version of the low-temperature rate factors R(s): 



where AV shear (s) is the activation volume required to nucleate a shear transformation. Our 
idea here was that, at temperatures well below the glass temperature, the transitions between 
STZ states are not thermally activated but, rather, are controlled entropically. That is, the 
rate factors are determined by the number of paths that the molecules within a zone can 
follow in moving around each other while going from one state to the other. The exponential 
factor in Eq. ()5.1|) is an approximation for a weighted measure of that number of paths. Its s- 
dependence means that greater weight must be given to paths moving in the direction of the 
stress than opposite to it. The exponential form of AV shear (s) is the simplest non-negative 
function that becomes arbitrarily small at large s and introduces just one new parameter, 
the effective STZ stiffness ft. The quasilinear version of the theory corresponds (roughly) to 
the limit of small s and small values of AVJ 3 s/iear /f/. 

Comparison of Eq. ([5.1)1 with Eq. ([3.1|) indicates that the natural way to include ther- 
mal effects in R(s) is simply to let Vf have the T-dependent Cohen-Grest form shown in 
Eq. (JH.2|) . This means that, at low T, the ratio AVQ hear /vf(T) becomes very large, which, 
in turn, implies that the functions C(s) and T(s) introduced in Eq. ([2.4)1 become strongly 
stress dependent, and the quasilinear approximations made in Eqs. ([2.8jl are no longer valid. 
Importantly, C(s) becomes very small for small s, so that plastic deformation is strongly 
suppressed at stresses appreciably below the yield stress. 

The strong stress dependence of C(s) and T{s) should be especially apparent in transient 
behavior of the kind shown in Figures 1, 2 and 7. Here, the initial response to loading at 
small stress will be almost entirely elastic, and plastic deformation will begin only later in 
the process. We shall have to use the fully nonlinear theory when undertaking more detailed 




(5.1) 
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comparisons with these kinds of experimental results. 

Shear localization: All of the analysis in this paper pertains to spatially homogeneous sys- 
tems. In order to make closer contact with experiments, we shall have to understand why 
and when these systems become unstable against shear banding and inhomogeneous failure 
modes, especially fracture. 

One mechanism for shear localization that we have not mentioned in this presentation is 
the elastic interaction between STZ's studied in As shown in that paper, an STZ-like 
event generates a quadrupolar stress field that induces other nearby events along preferred 
spatial directions and suppresses events elsewhere. The result is a tendency toward shear 
localization that should be interesting to examine in the context of this more general version 
of the STZ theory. 

A second mechanism that seems likely to play a role in shear localization is already built 
into our equations of motion when we write them in terms of spatially varying fields. From 
Eqs. (l3.7j) and (j3.10|) . we see that the STZ density A grows most rapidly, within limits, in 
regions where A already is large. This feedback effect, perhaps coupled to the effect of 
elastic interactions mentioned above, is our best guess at present about how shear banding 
will emerge in the STZ theory. 

Effective temperature and the interpretation of n eq : Finally, we return to Spaepen's sugges- 
tion that the density of STZ's might be directly related to the free volume as in Eq. (jl.3|) . 
Note that Vf appears in that equation, not as an extensive quantity related directly to the 
difference between the actual volume and some hypothetical close-packed volume, but rather 
as an intensive variable analogous to a temperature. For example, Mehta and Edwards}^ 
have introduced an intensive variable, thermodynamically conjugate to volume, which must 
govern density fluctuations in relations such as Eq. (jl.3|) or (|3.1|) in much the same way as 
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2911 has shown in a 



ordinary temperature governs energy fluctuations. One of us (MLF) 
recent molecular-dynamics simulation that while, as expected, the interior of a shear band 
has a high density of active STZ's, the actual free volume in that region is only slightly 
greater than in the rest of the system, but the intensity of density fluctuations - both 
compressive and dilational - is substantially higher. 

The idea that x = V f/V* m Eq. ()1.3|) might more generally be interpreted as an effective 
temperature seems especially appealing in light of our argument that the limiting steady- 
state value of A (or equivalently or n eq ) should depend upon the order in which we 
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take the limits T — > and e pl — > 0. More realistically, under steady-state conditions, x 
might approach some nonzero limiting value and remain there for indefinitely long times 
at sufficiently small T and for arbitrarily small but nonzero e pl . Conversely, x ( m suitable 
units) might approach the true temperature at large enough T or small enough e pl . 

There is increasin g ev idence that something like this happens in sheared foams or granular 
materials. (30, 31, 32, 3s| In those systems, the usual kinetic temperature is zero because the 
constituents have very large masses, but an effective temperature determined by response- 
fluctuation relations goes to a nonzero limit when the deformation rate becomes arbitrarily 
small. In our present system, there is a true kinetic temperature, but below the glass 
transition that temperature is so small that thermally assisted molecular rearrangements are 
effectively frozen out. During irreversible processes such as plastic deformation, therefore, 
the slow, configurational degrees of freedom characterized by \ might fall out of equilibrium 
with the fast, thermal (vibrational) degrees of freedom, and each may accurately be described 
by its own "temperature." We suspect that some such description of our system will be 
necessary in order to resolve the two-limit problem that we have encountered here. 



Acknowledgments 

J. Langer and L. Pechenik were supported primarily by U.S. Department of Energy Grant 
No. DE-FG03-99ER45762, and in part by the MRSEC Program of the National Science 
Foundation under Award No. DMR96-32716. M. Falk was supported by the National 
Science Foundation under Award No. DMR-0135009, and in part by the Dow Corning 
Foundation. 



[1] J.S. Langer and L. Pechenik, Phys. Rev. E 68, 061507 (2003). 
[2] D. Turnbull and M. Cohen, J. Chem. Phys. 52, 3038 (1970). 
[3] F.Spaepen, Acta Metall. 25(4), 407 (1977). 

[4] F. Spaepen and A. Taub, in Physics of Defects, edited by R. Balian and M. Kleman, 1981 Les 

Houches Lectures (North Holland, Amsterdam, 1981), p. 133. 
[5] A.S. Argon and L.T. Shi, Acta Metall. 31, 499 (1983). 

[6] D. Srolovitz, K. Maeda, V. Vitek and T.Egami, Phil. Mag. A 44, 847 (1981). 

26 



[7] D. Deng, A.S. Argon, and S. Yip, Phil. Trans. Roy. Soc. Lond. A 329, 549 (1989). 

[8] A.S. Krausz and H. Eyring, Deformation Kinetics, (John Wiley and Sons, New York, 1975). 

[9] P. De Hey, J. Sietsma, and A. Van Den Beukel, Acta Mater. 46, 5873 (1998). 

[10] M.L. Falk and J.S. Langer, Phys. Rev. E 57, 7192 (1998). 

[11] W.L. Johnson, J. Lu and M.D. Demetriou, Intermetallics 10(11-12), 1039 (2002). 

[12] H. Kato, Y. Kawamura, A. Inoue, and H.S. Chen, Applied Phys. Lett. 73, 3665 (1998). 

[13] J. Lu, G. Ravichandran, and W. L. Johnson, Acta Mater. 51, 3429 (2003). 

[14] L.O. Eastgate, J.S. Langer, and L. Pechenik, Phys. Rev. Lett. 90, 045506 (2003). 

[15] M.L. Falk, Ph.D. Thesis, Phys. Dept., Univ. of Calif., Santa Barbara (1998). 

[16] L. Pechenik, |cond-mat/ 0305516 

[17] J.S. Langer, Phys. Rev. E 64, 011504 (2001). 

[18] M.L. Falk and J.S. Langer, M.R.S. Bulletin 25, 40 (2000) 

[19] M.H. Cohen and G.S. Grest, Phys. Rev. B 20(3), 1077 (1979). 

[20] A. Masuhr, T.A. Waniuk, R. Busch, and W.L. Johnson, Phys. Rev. Lett. 82 2290 (1999). 

[21] A. Taub and F. Spaepen, Acta Metall. 28, 1781 (1980). 

[22] A. Greer and F. Spaepen, Ann. N.Y. Acad. Sci. 371, 218 (1981). 

[23] S. Tsao and F. Spaepen, Acta Metall. 33, 881 (1985). 

[24] PA. Duine, J. Sietsma and A. van den Beukel, Acta Metall. Mater. 40,743 (1992). 

[25] O.A. Hasan and M.C. Boyce, Polymer 34, 5085 (1993). 

[26] O.A. Hasan and M.C. Boyce, Polymer Eng. and Sci. 35, 331 (1995). 

[27] A. Mehta and S.F. Edwards, Physica A 157, 1091 (1990). 

[28] M.L. Falk, Y. Shi, "Strain Localization in a Molecular-Dynamics Model of a Metallic Glass," 
In " Supercooled Liquids, Glass Transition, and Bulk Metallic Glasses," Ed. T. Egami, A.L. 
Greer, A. Inoue and S. Ranganathan, (Mat. Res. Soc. Proc. Vol. 754, Pittsburgh, PA, 2003), 
CC6.20.1-6. 

[29] Y. Shi and M.L. Falk, unpublished. 

[30] I.K. Ono, C.S. O'Hern, D.J. Durian, S.A. Langer, A.J. Liu and S.R. Nagel, Phys. Rev. Lett. 
89, 095703 (2002). 

[31] L.F. Cugliandolo, J. Kurchan and L. Peliti, Phys. Rev. E 55, 3898 (1997). 

[32] P. Sollich, F. Lequeux, P. Hebraud and M.E. Cates, Phys. Rev. Lett. 78, 2020 (1997). 

[33] L. Berthier and J.-L. Barrat, Phys. Rev. Lett. 89, 095702 (2002). 

27 



